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WHAT IS CLAIMED IS: 




1 . A method for the computer animation of two colliding bodieg^y the 
iterative calculation of body positions from previous positions, said memod comprising: 

receiving data defining physical properties of saidjx5dies, said physical 
properties including positions and shapes of said bodies; 

determining a collision between said bocH 

calculating a closest-points vector bejween said two bodies from said 
physical properties of said bodies; 

calculating a contact force on^t least one of said bodies, said contact force 
along said closest-points vector and havijrfg a nonlinear relationship with respect to said 
closest-points vector such that said force increases sufficiently rapidly as said closest-points 
vector goes to zero to overcome inotions causing said collision between said bodies; 

calculating a position of said at least one of said bodies at a subsequent time 
interval in response to said/contact force upon said body; 

iteratively repeating said collision determining, closest-points vector 
calculating, contact ^orce calculating and position calculating steps for subsequent time 
intervals; and 

displaying calculated positions of said bodies at selected time intervals for a 
realistic animation of said colliding bodies. 



2. The method of claim 1 wherein said contact force F hasJJi©^elationship: 
F = ae b(c " x) -d, 

where a, b, c, and d are constants, x is tlje-iefigth of the closest-points vector. 



3. The mgthtfd of claim 1 wherein said contact force F has the relationship: 

F =a/(b-x) 

whepeif, and b are constants, and x is the length of the closest-points vector. 



4. The method of claim 1 wherein said bodies are defined with incompressible 
portions, said closest-points vector calculated between said incompressible portions of said 
two bodies. 



m 

29 



z> 

compi 



5. The method of claim 4 wherein said bodies are defined with compressible 
portions covering said incompressible portions, said collision determining step comprising 
determining whether said compressible portions of said bodies are in contact. 

P 

v 6. The method of claim 1 wherein at least one of said collidins / 6odies 



comprises a first body part connected to a second body part by a joinj^said joint having at 
least one rotational degree of freedom, a joint limit stop and a m^imum allowed rotation 
limit; and comprising: 

determining a rotation of said first body jpfot about said joint to reach said 
joint limit stop; 

calculating a torque upon said fij^t body part along said one rotational 
degree of freedom, said torque having a nonlinear relationship with respect to an angle of 
rotation from said maximum allowed rmation limit such that said torque increases 
sufficiently rapidly as said angle goes to zero to overcome motions causing said rotation of 
said first body part about said joint; 

calculating a potational position of said body part at a subsequent time 
interval in response to saja torque upon said body part; and 

iterativ^ly repeating said rotation determining, said torque calculating, and 
said rotational position calculating steps for subsequent time intervals. 

7. The method of claim 6 wherein said torque T has the relationship: 

T = ae b(c - e) -d, 

where a, b, c, and d are constants, and 6 is an angle of said body part from said maximum 
allowed rotation limit. 

8. The method of claim 6 wherein said torque T has the relationship: 

T =a/(b-0) 

where a and b are constants, 0 is an angle of said body part from said maximum allowed 
rotation limit. 
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1 9. The method of claim 1 further comprising: 

2 determining a fault in any one of said closest-points vector, contact force, 

3 and position calculating steps; 

4 reinitializing said faulty calculating step; and 

5 performing said reinitialized calculating step with a smaller time interval. 

1 10. The method of claim 9 further comprising: 

2 repeating said fault determining, reinitializing and reinitialized calculating 

3 step performing steps at smaller and smaller time intervals until said fault is eliminated. 



1 11. The method of claim 9 wherein said smaller time interval is one-half of a 

2 2 previous time interval. 

*F 1 12. The method of claim 9 wherein said fault comprises floating-point numeric 

S 

j: 2 overflow. 

£3 1 13. The method of claim 9 wherein said fault comprises said closest-points 

p 2 vector being less than zero. 



1 14. The method of claim 9 wherein said fault comprises physically incorrect 

2 collision situations involving excessive interpenetration of bodies. 

1 15. The method of claim 6 further comprising: 

2 determining a fault in any one of said closest-points vector, contact force, 

3 position, torque and rotational position calculating steps; 

4 reinitializing said faulty calculating step; and 

5 performing said reinitialized calculating step with a smaller time interval. 

1 16. The method of claim 15 wherein said joint comprises a gimbal joint and said 

2 fault comprises said gimbal joint reaching a gimbal lock singularity. 
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17. The method of claim 15 wherein said rotational position calculating step 
comprises calculating Euler parameters and said fault comprises non-normalized Euler 
parameters. 




18 The method of claim 15 wherein said fault comprises first body part rotates 
beyond said maximum allowed rotation limit. 

19. A computer program for the realistic computer animation of twef colliding 
bodies by the iterative calculation of body positions from previous positimr^, comprising: 

code that receives data defining physical properties ofjmd bodies, said 
physical properties including positions and shapes of said bodies^ 
code that determines a collision between said Wxlies; 
code that calculates a closest-points vectoj/tfetween said two bodies from 
said physical properties of said bodies; 

code that calculates a contact force/6n at least one of said bodies, said 
contact force along said closest-points vector amd having a nonlinear relationship with 
respect to said closest-points vector such tjzat said force increases sufficiently rapidly as 
said closest-points vector goes to zero yf overcome motions causing said collision between 
said bodies; 

code that calculates/a position of said at least one of said bodies at a 
subsequent time interval in response to said contact force upon said body; 

code that iteraively repeats said collision determining, closest-points vector 
calculating, contact forc^/calculating and position calculating steps for subsequent time 
intervals; 

code ihat allows a computer displa y to disp lay„calculated positions of said 
bodies at selected /ime intervals for a realistic animation of said colliding bodies; and 
computer readable medium that stores said codes. 



20. The computer program of claim 19 wherein said contact fore 
relationship: 

F = 

where a, b, j^amTcfare constants, x is the length of the closest-points vector. 
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1 22. The computer program of claim 19 wherein said bodies are defined with 

2 incompressible portions, said closest-points vector calculated between said incompressible 

3 portions of said two bodies. 



1 23. The computer program of claim 22 wherein said bodies are defined with 

2 compressible portions covering said incompressible portions, said collision determining 

3 step comprises code that determines whether said compressible portions of said bodies are 

4 in contact. 



24. The computer program of claim 19 wherein at least one of said colliding 
todies is defined as a first body part connected to a second body part by a joint, said joint 

-f" 3 having at least one rotational degree of freedom, a joint limit stop and a maximum allowed 

4 rotation limit; and said computer program comprising: 

5 code that determines a rotation of said $frst body part about said joint to 

6 reach said joint limit stop; 

7 code that calculates a torque upon s&id first body part along said one 

8 rotational degree of freedom, said torque having a nonlinear relationship with respect to an 

9 angle of rotation from said maximum allowedr rotation limit such that said torque increases 

1 0 sufficiently rapidly as said angle goes to z^ro to overcome motions causing said rotation of 

1 1 said first body part about said joint; 

12 code that calculates a rotational position of said body part at a subsequent 

13 time interval in response to said torque upon said body part; and 

14 code that iteratively repeats said rotation determining, said torque 

15 calculating, and said rotational Position calculating steps for subsequent time intervals. 
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25. The computer program of claim 24 wherein said torque T has the 
relationship: 

T = ae b(c - e) -d, 

where a, b, c, and d are constants, and 0 is an angle of said body part from said maximum 
allowed rotation limit. 

26. The method of claim 24 wherein said torque T has the relationship: 

T =a/(b-0) 

where a and b are constants, 0 is an angle of said body part from said maximum allowed 
rotation limit. 

27. The computer program of claim 19 further comprising: 

code that determines a fault in an operation of any one of said closest-points 
vector, contact force, and position calculating codes; 

code that reinitializes said faulty operation; and 

code that performs said reinitialized operation with a smaller time interval. 

28. The computer program of claim 27 further comprising: 

code that repeats operations of said fault determining, reinitializing and 
reinitialized operation performing codes at smaller and smaller time intervals until said 
fault is eliminated. 

29. The computer program of claim 27 wherein said smaller time interval is 
one-half of a previous time interval. 

30. The computer program of claim 27 wherein said fault comprises floating- 
point numeric overflow. 

31. The computer program of claim 27 wherein said fault comprises said 
closest-points vector being less than zero. 
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1 32. The computer program of claim 27 wherein said fault comprises physically 

2 incorrect collision situations involving excessive interpenetration of bodies. 

1 33. The computer program of claim 24 further comprising: 

2 code that determines a fault in an operation of any one of said closest-points 

3 vector, contact force, position, torque and rotational position calculating codes; 

4 code that reinitializes said faulty operation; and 

5 code that performs said reinitialized operation with a smaller time interval. 
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1 34. The computer program of claim 33 wherein said joint comprises a gimbal 

2 joint and said fault comprises said gimbal joint reaching a gimbal lock singularity. 

1 35. The computer program of claim 33 wherein operation of said rotational 

2 position calculating code comprises calculating Euler parameters and said fault comprises 

3 non-normalized Euler parameters. 
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part rotating beyond 



The computer program ofclaiiii^t5--^Kerein said fault comprises first body 

allowed rotation limit. 
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